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NATION ^ The in a previous paper for extending the neutron scattering formalism of Zemach 

I • l *~ "and Glauber to any type of molecular rotator have been employed to derive generalized forms of the 

SPACE adm: dh^^ep^TiQN sections for rotator scattering. A mass- ratio expansion for the treatment of the high-energy 
“ 1 ' limit is illustrated on the classical cross section and then employed in the treatment of the more general 

Lancl-y fiL quAhtditf-mechanical expression for the differential cross section. The results apply to an arbitrarily 
asymmetric rotator. The very low energy approximation is carried out for the symmetric rotator, and 
the procedure is compared with the explicit summing of the partial cross sections for individual rotational 
transitions. The inelastic correction to the static approximation for interference scattering is calculated 
to an accuracy of first order in the mass ratios for the case of the symmetric rotator. 
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1. INTRODUCTION 

T HE formalism of Zemach and Glauber 1 for 
treating the scattering of low-energy neutrons 
by chemically bound nuclei in a molecule was applied 
by them only to the simplest types of molecules as 
far as the effects of the rotational degrees of freedom 
are concerned. Procedures for generalizing their 
results to any type of molecular rotator have been 
given recently 2 and were used to calculate the distribu- 
tion of scattered neutrons over energy and angle up 
to the accuracy of the first quantum-mechanical 
correction to the classical cross section. The remaining 
generalized differential cross sections for rotator 
scattering corresponding to the other types of approxi- 
mations originally examined by them have now been 
derived and will be presented here. 

In the Fermi pseudo-potential approximation the 
differential cross section of the entire rotator for a 
given initial rotational state can be written in the 
notation of reference 2 as 

*«(#,€)- (*/2rto)£ f (1.1) 

•'-oo 

in which the expectation value in the integrand has 
the form 

<*,',)= exp(iic-b,0<r tin exp(-iVb,) | ft). 

Into this basic expression derived by Zemach and 
Glauber there has been incorporated by the use of 
operator techniques the summation over all final 
states allowed by the conservation laws. The total 
potential presented to the neutron by the rotator is 
composed additively of the individual nuclear potentials, 
so that Eq. (1.1), which is based on the first Born approx- 
imation, has the form of a sum of terms which are of two 
general types. A direct scattering term is one which is 

‘A. C. Zemach and R. J. Glauber, Phys. Rev. 101, 118, 129 
(1956). 

* H. C. Volkin, Phys. Rev. 113, 866 (1959). This paper will 
hereinafter be referred to as I and an equation appearing in it 
will be designated by the numeral I preceding the equation 
number. 


attributable to a single nucleus of the rotator while ai 
interference term involves two scattering nuclei. The dif 
ferential cross section for the direct scattering from a nu 
cleus in the classical or high-energy limit is developed in 
powers of the ratio of the neutron mass to an effective “ro- 
tational mass” of the molecule. The calculation is carried 
through up to terms of the fourth order in the mass 
ratio. The expansion of the general expression for the 
cross section in powers of mass ratios is then discussed 
and the calculations up to terms of the second order are 
presented. The very low energy region is considered 
next and an approximate cross section which is suitable 
for this limit is derived for the case of the symmetric 
rotator. When the neutron energy is very low, the 
number of energetically allowed rotational transitions 
may be so limited that a direct summing of the partial 
cross sections for individual transitions is feasible. 
This alternative procedure is compared with the 
preceding calculation where the low-energy approxima- 
tion was applied to a direct scattering term of Eq. (1.1), 
which contains implicitly the sum over all possible 
transitions. Finally the portion of the differential cross 
section which arises from the interference effects due 
to the presence of more than one scattering nucleus in 
the rotator and for which no classical approximation 
exists is treated by means of the mass ratio expansion 
and the calculation for the symmetric rotator is 
exhibited up to the first order terms. 

2. RESUME OF BASIC FORMULAS 

The notation to be used here is taken over directly 
from I. We consider the scattering of a neutron with 
initial momentum k 0 to final momentum k by a rotator 
with fixed center which is in the initial state \JM) 
having energy Ej , angular momentum /, and compo- 
nent of angular momentum along an axis fixed in 
space M. The third quantum number required to 
specify a rotator state is omitted for simplicity in the 
case of an asymmetric rotator energy eigenstate, where 
it has no relevant physical significance. Scattering with 
momentum gain K=k— k 0 to the neutron corresponds 
to a certain scattering angle# and energy gain e= (2 m) -1 
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X(k 2 —k 0 2 ). The general expression for the differential 
cross section of the entire rotator is given by Eq. (1.1). 
The direct scattering term due to the vth nucleus is 
the term with /= v , while the two terms involving v f 
and v, with y V v, are the interference terms contributed 
by the corresponding pair of nuclei. The states which 
differ only in orientation of the angular momentum 
vector, i.e,, states of different M, are of course degen- 
erate in energy. The cross sections must therefore be 
averaged over the initial values of M. Also it is sufficient 
to consider the scattering as spin independent, since the 
consequences of spin dependence have been treated 
explicitly by Zemach and Glauber. Then the part of 
the differential cross section which is due to a single 
nucleus having position vector b relative to the center 
and scattering length a is ^ 

<rj(#,€) = k(27rk 0 )- 1 f dte~ i(i (x)j , (2.1) 

— 00 

where the expectation value is given by 


for which M specifies the component of angular momen- 
tum. Then a complete set of commuting observables 
for the rotator is L 2 , the component of L along a body 
axis, say 7 r , and the component of L along a space 
axis, say L z . The corresponding set of quantum 

numbers 7, K , and M label a complete orthogonal set 
of basic states | JKM) for the system. An energy 
eigenstate | JM) can always be expressed as 

\JM)~ £ a K \JKM). (2.4) 

K=—J 

It is to be understood that in the expansion the compo- 
nent K refers to the z axis of the body system chosen 
to resolve components of the vector and tensor quanti- 
ties. Then the operation of L in (2.2) is specified by 

L Z \JKM)=K\JKM\ (2.5a) 

(L X ^L V ) | JKM)= [(7=Ftf ) (J±K+ 1)]» (2.5b) 

X| 7, K±1,M). 


<x)j=a 2 (2/+l)'~ 1 £ (JM\e-^ A + B +^\JM). (2.2) 

M=—J 

The Hermitian operators A, B,C are defined by 
A=H-Ej, 

B= — k*Q* L+^ic*K 
= — L-Qt-K— §fC-K, 

C=§K‘R-k. 

The rotator Hamiltonian may be written as H= |L • G 
*L, where L is the angular momentum operator. The 
tensor G is the inverse of the moment of inertia tensor, 
G=I“ 1 . We further define the tensor quantities By 
= timjbm, R=BtGB, Q=BtG, and the vector c m =biGi m 
— (Trace G)6 m . 

In evaluating the matrix elements (2.2) it will be 
necessary to express the vector and tensor quantities in 
terms of components relative to some system of body 
axes fixed in the rotator, but not necessarily the system 
of principal axes. When so taken the components of 
the vector b and of the tensors given above are constants 
and commute with the components of L and k as well 
as with each other. The momentum transfer vector k, 
which is fixed with respect to space axes, becomes a 
dynamical variable when seen from the rotating system. 
The components of k and L may easily be shown to 
obey the following commutation relations in this 
description. 

ZLM = iUikLt, — ieykKic, 

[W>0. (2.3) 

Let us now designate the set of body axes along which 
components are being taken by xyz and a set of space 
axes by XYZ , one of the latter, say Z, being taken for 
convenience as parallel to k and as the space direction 


In terms of the basic set of states we may write the 
expectation value (2.2) in the form 


with 


(x)j=a 2 £ <xk*&k'{x)jkk'> 

K,K' 


(x)jkk’= (27+ 1)- 1 Em (JKM | £-«<a+b+o | JK'M). 

( 2 . 6 ) 


The Euler angles which specify the orientation of the 
body axes relative to the space axes are defined as 
follows: 6 and <p are the polar and azimuthal angles, 
respectively, of the z axis and \p is the spin angle about 
this axis. The wave function representatives of the 
basic set in these variables may be specified by 

= [(27+ 1)/8 

where /x = c os0 and is the KMt\i matrix element 

of the 7th irreducible unitary representation of the 
rotation group. These matrix elements of the rotation 
operator are given explicitly as 

UM tm =(jk \uwm\jm), 

with 

V (vfljfr) = e ^ Js e iej2 e i<pJi , 

where the 7 t satisfy the usual angular momentum 
commutation relations and | jm) is an eigenstate of J 2 
and Jz. Because of the different conventions which 
have been used by various authors we mention that the 
rotation operator U given here transforms an arbitrary 
state vector (£'| and arbitrary operator O in the 
following manner: (£' | U (a,P,y) and U~ l OU are the 
state vector_and operator which look in the reference 
system XYZ produced by the rotation ( aj3y ) from the 
original reference system XYZ as (£' | and 0 do, respec- 
tively, in XYZ. In terms of the spherical harmonics 
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y im defined by Condon and Shortley 8 it may be shown 
that 

U^mWm = 0/(2/+ l)]»F«(fl ><9 ), 
U^ko(vM)= (-l)*0/(2/+l)]*F,*(W). (2 ' 7) 

For certain purposes it will be advantageous to use the 
principal body axes xyz in which case the tensor I is 
diagonal. We write the expansion (2.4) explicitly in 
this case as 

\JM)= Y,R&r\JKM), (2.8) 

where K is the component of L along the principal 
z axis. The form (2.8) is the one generally presented in 
the literature. It was shown in Appendix C of I that 
the expansion coefficients are related by the transforma- 
tion law 

<XK = JllcaRU( J) KK (a,P,y ) , (2.9) 

where a, /?, y are the Euler angles of the principal axes 
system relative to the system xyz. The principal axes 
have the advantage that for a symmetric rotator the 
states | JKM) are themselves the energy eigenstates if 
the z axis is the axis of symmetry, i.e., under these 
conditions _the _expansion (2.8) reduces to az=l, 
aR, = 0foTK'^K. 

3. CROSS SECTIONS FOR THE DIRECT 
SCATTERING 

For large values of k 0 b the term C in the exponential 
operator of (2.2) dominates. The terms A and B are 
of order (i <b)~ 2 and ( Kb)~~ l relative to C. If we write 
\=—it and the exponential operator as exp(XC)/(X) 
with /(X) = exp(-XC) expX(^+-B+C), the classical 
approximation corresponds to replacing /(X) by the 
unit operator, while the first quantum-mechanical 
correction consists of retaining terms in the expansion 
of fOO U P to order (k£)“ 2 . These cross sections are 
derived in I, Sec. 3. The result of the classical approxi- 
mation can be written as 

<Tci (#, c) = 2a 2 k\jr (rir 2 ) £x~y i K (p) 

for — £r 2 /c 2 <€< 0, 

= 2 a?k[jr (r ir 2 ) *k#c t x~y i qK (q) (3.1) 

for — ^riK 2 < € < — £r 2 x 2 , 

= 0 otherwise, 

where r\ and r 2 (ri>r 2 >0) are the two nonvanishing 
eigenvalues of R, K is the complete elliptic integral of 
the first kind, and 

as 2 — l + 2e(r 1 K 2 ) -1 , 
f=q~ 2 =r~ 1 (cr 2 ~~ 1 ), 
r=r 2 (ri-r 2 )- 1 . 

The values of rj and r 2 are given explicitly in terms of 
the principal moments of inertia and the components of 
b relative to the principal axes in I (3.2). 

8 E. U. Condon and G. H. Shortley, Theory of Atomic Spectra 
(Cambridge University Press, Cambridge, 1935), p. 52. 


On introducing the dimensionless variables 
A = ^^ 0 “ 1 Cl+5(l — kokr 1 cos??)] -1 , 

5 = tut i ( 1 — ir 2 ) , 

we find the differential cross section a c i(t?) = J ' (r(&,e)de 
becomes 

<r c i(t?) = 2a 2 (r 1 r 2 ) i [ir(r 1 — rj)]- 1 ^ dp pK(p ) 

XMl+rf)~i+J^ <fy ? JT( ? )A(r+tf)-»J. (3.2) 

With the help of the relation M 0 ' l = (1+5)— ^5 cost? 
+ (1— 5 2 sin 2 t?)*], we may express A in terms of the 
variable 5 alone and can then carry out the power 
series expansion 

A= 2Z d n $ n . 

n=0 

The first few expansion coefficients are 

da— 1, di= — 2(1 — cos??), 

<* 2 =i(l-cosj?)( 5-3 cost?), ^ 3 = — 3 (1 — cost?) 2 , 
d 4 =|(l — cost?) 2 (27— 10 cost?— 5 cos 2 t?). 

As a function of p the variable 8 has the form 8 = mrirp- 
X(1 +rf)~ l while as a function of q it takes the form 
8=mrir(r-\-q 2 )~ l . 

If we insert the series expansion of A and use the 
appropriate functional form of 8 in the two integrals of 
Eq. (3.2), we are led to the relation 

<Tci(t?)= (2a 2 / 7r)[r(l+r)]l d n (—2mrir) n 

XHI-'-C 2»+l)- l ]0/ar)»;(r), 

with g(r) = h(r)+r-ih(r~ l ) where h(r) = fo 1 dq qK(q) 
X(r+? 2 )H. 

In order to evaluate the integral h(r) one can derive 
from certain differential relations for the complete 
elliptic integrals the indefinite integral 4 

rj' dq qK(q)(r+q 2 )~f=q 2 K(q)(r+q 2 )~i 

-n(r; 9 )(r+f)*, 

where II is the complete elliptic integral of the third 
kind. Substituting the desired limits yields 

h(r) = 0(l+r)]-i[i,r— tan" 1 (»•♦)], 

so that 

gW = §»[f(l+f)]-*. 

We obtain finally the expansion of the classical cross 

4 Integrals similar to this have been calculated in full detail by 
K. F. Miiller, Arch. Electrotech. 17, 336 (1926). 
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section in powers of the mass ratio mri : 

o' ci (#)/a 2 =1+| (wr 0 ( 1 + 2r ) ( 1 + f ) _1 #i+ (1/15) (mr i) 2 
X (3+8r+8r 2 )(l+r)~ 2 ^2+ (1/35) (wfi) 3 
X (5+18f+24r 2 +16r 3 )(l+r)-^ 3 
+ (1/315) (wfi) 4 (35+ 160r+288r 2 

+ 256f»+128f 4 )(l+r)- 4 d 4 +- • ■. (3.3) 

The mass ratio expansion for the general expression 
(2.2) will now be considered. Following Zemach and 
Glauber we take the Taylor expansion of the exponential 
operator and obtain thereby a corresponding expansion 
of the expectation value in powers of (if) which we 
write as 

(x)j= £ «»(*'<)"• 

n=0 

This is substituted into (2.1) and then the integration 
over the variable e is performed. The differential cross 
section takes the form of a series expansion in powers of 
the mass ratios (mb x 2 /I >), 

aj(d) = k o -1 En [(mk-'d/dkyken'lk^ko. (3.4) 

In the expansion of expQ— it(A-\-B+C)~\ we see 
immediately that terms having A all the way to the 
left or right give no contribution to the expectation 
values in Eq. (2.2) since A\JM)=0. Terms of odd 
degree in B likewise contribute nothing to the expecta- 
tion values. It may be seen from Eq. (2.6) that they 
lead to integrals of the form (.44) which vanish since 
their integrands are odd. When we retain only contribut- 
ing terms the expansion of (2.6) to third order in t has 
the form 

(x)jKK' = bKK'A- (2/+1) -1 YLm (JKM\ (— . it)C 

+ ±(ity(B*+c>)-(i/6)(ity 
X (BAB+CAC+B 2 C+CB 2 +BCB+C*) 

X| JK'M). (3.5) 

One choice of body axes that suggests itself for 
evaluating the matrix elements (3.5) is the system in 
which R is diagonal. The position vector b is an 
eigenvector of R corresponding to the eigenvalue zero, 
and we take the direction of the z axis along this vector. 
The components of 0, G, and c in this system are given 
in terms of r h r 2 and the principal moments of inertia 
in I, Sec. 3. 

The procedure for evaluating the expectation values 
(2.6) is first to commute the angular momentum 
components all the way to the right or left of the terms 
in which they appear. Their operation on the basic 
states is described by Eqs. (2.5). We then go over to 
the Euler angle representation. Since no factors depend- 
ent on M have been introduced, the average over M 
can be carried out immediately by using the unitary 
property of the wave functions. This is illustrated in 
the Appendix where the resulting integrals are eval- 
uated. Finally the expansion coefficients e n are inserted 


into Eq. (3.4). We note that the contributions of C 
and C 2 in Eq. (3.5) are formally independent of the 
choice of body axes. This can be seen as follows. Since 
C does not contain the angular momentum operator, 
an arbitrary function f(C) does not connect states of 
different K regardless of the axes system in which 
components are taken. More precisely we have that 

(2/+1)- 1 Em (JK'M\f(C)\JKM) 

= (1/8 f{C)du, (3.6) 

with the integral on the right-hand side being taken 
over all orientations of the rotator. Since C is a scalar 
invariant, its value for a given orientation of the 
rotator is independent of the choice of body axes, and 
therefore the integral is likewise. Using axes which 
diagonalize R and the evaluation of integrals in the 
Appendix, we easily find that the right-hand side of 
Eq. (3.6) for the case f(C) — C equals i+f 2 ) 

and for f(C) = C 2 equals (l/60)5iCK'/c 4 (3ri 2 +2rir 2 +3r 2 2 ). 
The complete computation to second order in the mass 
ratios yields 

a j (#) /a 2 — 1 — f m (r 1 + r 2 ) ( 1 — cos#) + ( 1/30) m 2 

X (3ri 2 +2fir 2 +3r 2 2 )(l — cos#) (5 —3 cos#) 
+|w 2 (W)- 2 {Kn 2 +r 2 2 )[/(/+l)- (K 2 )-] 
+Kn+r2) 2 +6 4 (G 13 2 +G 2 3 2 )[<i^)+i] 

+E* l(J-K+l)(J+Kmb 2 r 2 Gn(2K-\) 

X0t(«irWi)+i(r2 2 -ri 2 )(/-^+2)i. 

X (J-\~K— l)*(R(aK*atf_ 2 ) 

+ib 2 r 1 G^(2K-\)g(a^a K - l )']), (3.7) 

where (i£ 2 ) = Ex (*k*(xkK 2 and 6t and $ denote the 
real and imaginary parts, respectively. 

In the case of the symmetric rotator the symmetry 
about the principal z axis permits placing the nucleus in 
the plane of the principal x and z axes. The components 
of b relative to the principal axes are then (&i,0,& 3 ). Thus 
the z axis lies in the xz plane andmakes an angle j3o=tan~ 1 
y>(bi/bf) with the symmetry axis. The Euler angles of the 
principal axes relative to xyz are seen to be (0, — 0o, 0). 
For the energy eigenstate \JKM), we have in Eq. 
(2.8) that 

a} c=l, ctK f = 0 for K f 7*K. (3.8) 

The coefficients a K are then given by the transformation 
law (2.9) as 

OLK=U^fCK( 0, — 00,0). (3.9) 

Formally an average over K and —K is required, but 
the cross section for the symmetric rotator is in general 
even in K , since the states | JKM) and \J — K —M) 
have the same value for the matrix element in Eq. (2.2) 
and an average over all values of M is taken. That the 
expression (3.7) for <jj{0) is even in K may be verified 
explicitly by using the relation 

U^-r k ( 0, -/So, 0)= (~l) K + R U iJ) R- K (0, — /3o, 0). 
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In I the expression for the first quantum correction 
to the classical cross section was based upon the use of 
body axes which diagonalize R. It was erroneously 
asserted there (p. 871) that the special case of a sym- 
metric rotator corresponded to the condition (3.8) on 
the ax. The coefficients given by (3.9) are instead the 
generally valid ones. Only when the scattering nucleus 
lies on the symmetry axis is (3.8) also applicable to the 
ax, for then U {J) rk( 0, 0, 0) = 5^x. The results for the 
examples of symmetric rotators discussed in I remain 
unaffected, since for the spherical and linear rotators 
the axes xyz defined above are also principal axes. 

An alternative form for the mass ratio expansion of 
the cross section may be obtained by using the principal 
body axes xyz to evaluate the matrix elements (3.5). 
Again the results will be given to second order in the 
mass ratios. When the principal body axes are used, G 
is diagonal and in terms of the principal moments of 
inertia is given by Ga=Ic l - The quantities appearing 
in the operator B then take the form 


Cm b m (G, 

nm Trace G), 

(3.10) 

0 

bzlf 1 

-bilr 1 ) 


-Wf 1 

0 

bJr 1 

, (3.11) 

bJr 1 

-bih~' 

0 J 



where the components bi are, of course, taken with 
respect to the principal axes. The same procedure used 
to obtain (3.7) is employed to evaluate the contribution 
of B 2 in Eq. (3.5). 

The alternative expression which is found for the 
differential cross section by the means described is 
perhaps a more useful one than Eq. (3.7). The expansion 
coefficients olr are the ones generally evaluated in 
computations of molecular wave functions and they 
possess the more elementary form (3.8) when the 
molecules are symmetric. The differential cross section 
may be written as 

c Tj{d)/cP 

= \ — \m (Trace R)(l — cost?) + (1/30) w 2 
X (3ri 2 +2rir 2 +3r 2 2 )(l — cost?) (5—3 cost?) 

+ > 2 ^0- 2 (§[/r 2 (6 2 2 +i3 2 ) + /2- 2 (&X 2 +^3 2 )] 

X[/(/+ 1)- (K*)-]+I z ->(bS+bS)(K>) 

+llbxKh- 1 +h- 1 ) 2 +b2 2 (ir 1 +h- 1 ) 2 

+W(/r 1 +/2- 1 ) 2 ]+L^[(/-X+i)(7+X)]i 

* X { - (2 K - 1 j/r’iaC/f (R(a**ajn) 

+t/ 2 - 1 ^^(a^*a/?-i)]+§[/i- 2 (62 2 +&3 2 ) 
-/2“ 2 (6i 2 W)][(/-X+2)(/+K-l)]i 
X (R (a£*a£_ 2 ) — il f" 1 /^ i —1 4 1^2 
X[(/-X+2)(/+Z-l)]^(a^_ 2 )}). (3.12) 

In order to specialize (3.12) to the case of a symmetric 
rotator, let us take the z axis as the axis of symmetry. 
Then /i— / 2 and the symmetry permits taking & 2 = 0. 
Furthermore the condition (3.8) is satisfied, so that the 


summation term vanishes and (K 2 )=K 2 . The resulting 
expression is even in K . With /i=/ 2 =/ and b 2 =bi 2 +b z 2 
the differential cross section for the symmetric rotator 
becomes 

<jj{d)/a 2 

= l-§ (■ m/IhWh+bSI+bNd (1 - cost?) 

+ (1/ 30) (mf 27 3 ) 2 [3& 4 / 3 2 + 2& 2 / 3 (b\ 2 I + b z 2 1 3 ) 

+ 3 (5 1 2 /+ 6 3 2 / 3 ) 2 ] ( 1 — cost?) (5—3 cost?) 

+IWW/ 3 ) 2 {|/ 3 2 (^i 3 2 )[/(/+l)-^] 

+Pb l 2 K 2 +ib l 2 {I+h) 2 +b z 2 I z 2 }, (3.13) 

Two special cases of symmetric rotators are of 
particular interest. For the spherical_ rotator we set 
I Z =I and average over all values of K, which has the 
effect of replacing K 2 by ^/(7+l), obtaining 

aj(#)/a 2 = 1 -%{mb 2 /I) (1 - cost?)+ (4/15) {mb 2 / 1) 2 
X (1 — cost?) (5—3 cost?)+^(fw6/^o/) 2 

X [§/(/+ 1)+1]. (3.13a) 

For the linear rotator we need only set K=0 and 
&i=0 in (3.13). The differential cross section for a 
linear rotator is thus 


aj (t?)/ a 2 =1-| {mP/I) (1 - cost?) + (4/15) (mb 2 / 1) 2 
X (1 — cost?) (5— 3 cost ))+%(fnb/koI) 2 

X[/(/+l)+l], (3.13b) 


which agrees with the results stated by Zemach and 
Glauber. 

We consider next the region of very low neutron 
energies where the preceding approximations are 
inadequate. Here the mass-ratio expansion may be 
expected either to converge slowly or, for neutron 
energies less than the rotational level spacing, to 
diverge. The low-energy approximation is based on 
the property of Eq. (2.2) that in this limit the term A 
dominates in the exponential operator. In the approxi- 
mation to order {kjb) 2 the operator is expanded and 
terms up to second order in B and first order in C are 
retained. The desired form of the operator can then be 
written as 


e X(X+B+C) =1+X ( 5+C ) + 



X' 


d\" B(T A B 


+A(---)+{---)A+0(BC). (3.14) 


As explained above the terms with A all the way to 
the left or right give vanishing contributions in (2.2) 
and the term linear in B does likewise in the matrix 
elements (2.6). The matrix elements of C between states 
| JKM) are known from the previous calculations. In 
order to handle the integral term explicitly, we now 
restrict our considerations to the case of a symmetric 
rotator. The symmetric rotator Hamiltonian and 
consequently the operator A are diagonal in the basic 
set of states for which the body-axis component of 
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angular momentum refers to the symmetry axis. 
Moreover the eigenvalues are known exactly. Calling 
the symmetry axis the z axis and the common moment 
of inertia about the other two principal axes /, we 
write the well-known result 

H\JKM)=E{JS)\JKM) 

E{jS) = hI~ l J{J+^)+Wz^-I~ l )K\ (3.15) 

To find the contribution of the integral term in 
(3.14) components will be taken in a system of principal 
axes bearing the notation just given. The symmetry 
about the z axis further allows us to select the xz plane 
so as to contain the scattering nucleus, and hence we 
may set bi — 0 in Eqs. (3.10) and (3.11). Since the 
body axes thus specified are the only ones entering the 
calculations, we shall for simplicity omit the bar on K 
which has heretofore been associated with the use of 
( 2 . 8 ). 

Let us specify the initial state of the rotator by 
\JKM). In order to evaluate the expectation value 
(JKM | B exp(\"A)B\JKM), we determine the opera- 
tion of B on the state vector. Once the vector B \ JKM) 
is known we need only operate on it with the exponential 
operator using (3.15) and then take the scalar product 
of the result with the adjoint vector since B is Hermi- 
tian. The desired evaluation can be effected in the 
explicit representation given by the Euler angles 
of the body axes. The components of k 
relative to these axes are given by 

k\ = — k sin0 cos^ 

= - (2)“*k[£/«> 10 (co)- <»_ 10 («)], 

K2=Ksindsin\p (3.16) 

= — ^(2)“ i /c[Z7 (l ) 10 (w)+ £/ (1) _io(co)], 
kz=k cos0=/c{/ (1> oo(oo). 

The operation of B on the symmetric top wave function 
N jU {J) km(u>) is found as follows. Insert (3.16) into one 
of the forms given for B . Then the products of matrix 
elements of U are reduced by means of the tabulated 
Clebsch-Gordan coefficients (JKlk\JljK+k) and the 
operation of the body-axis components of L is given by 
(2.5). The relevant reduction formula 5 is 

U (l) k o (co) ZJ( J) km (to) 
j+i 

= E (JKlk\JljK+k)U^K +ktM (u) 

• . MJ-l\ 

X(JMlO\JljM). (3.17) 

After performing the computations outlined and 
averaging over M, we collect the contributions of the 
various terms in (3.14) to the expectation value (2.2) 
and obtain finally for the scattering cross section of a 

6 See, for example, A. R. Edmonds, Angular Momentum in 
Quantum Mechanics (Princeton University Press, Princeton, 
1957), p. 60. 


symmetric rotator in this approximation 
(#,€)/ a 2 

= ttrtl- iW+WJ-'(J+ l)- l K*y(e) 
+Pfc<fW(2/+ 1) _1 { (/+ 1 - K) (. J+l+K ) 

X (/ +l)“ 1 5[e+A(/+l, K)~\ 

+ (J-K)(J+K)J-'8le+A(J-l, IE)]} 

+ 1) (2 /+ 1) J-i 

X{(J+l+K)(J+2+K)8le+A(J+l, K+ 1)] 

+ (/+1- j K)(/+2-JE)«[6+A(/+1 > K- 1)]} 

+ [/(/+ 1)]" 1 * (J-K) (7+ 1 +K ) 

XS[e+A (/, k+1)]+(J+K)(J+1-K) 

X6[e+A (/, X-l)]} + [/(2/+l)]-i 

X{(J-K)(J-l-K)8te+A(J- 1,JC+1)] 

+ (J+K ) (/- l+2C)fi[e+ A(/- 1, K- 1)]}), (3.18) 

where the rotational level spacing is written as A(J h Ki) 
= E(J h K 1 )—E(J,K). Writing p= - (2m/k 0 2 )A(J h Ki) 
we have for the transition (/,7T) — > (Ji,Ki) 

kk 0 ~ l K 2 = 2^ 0 2 (l+p) i [l+^p— (l+p)* cos??]. (3.19) 

The condition |p|<Kl holds when the mass ratios are 
sufficiently small. We can expand the right-hand side 
of Eq. (3.19) in powers of p and obtain thereby a 
mass ratio expansion of the low-energy cross section 
(3.18). Retaining terms to second order in p yields the 
first, second, and fourth of the four terms presented 
in the more general expansion (3.13). This agrees with 
the observation of Zemach and Glauber that the mass 
ratio expansion of the low-energy approximation 
contains the terms which dominate at low energies in 
the more general expansion (3.5). 

It is of interest to consider the explicit summing of 
the partial cross sections for individual rotational 
transitions. The procedures involved will be described 
for the symmetric rotator and the very low energy 
case will then be compared to the preceding calculation. 
We begin with the expression for the differential cross 
section of the scattering process, in which the rotator 
undergoes a transition from the initial state ypi to the 
final state ^7, 

<r f »(??) = E a v '(iy(k/ko)M V '*M, n (3.20) 

V* , V 

where the sum is taken over all nuclei of the rotator and 
M , = (if'/ 1 exp ( - Ik ■ b „) | \pi) . 

The complete cross section for a given initial state is 
the sum of the partial cross sections (3.20) over all 
final states allowed by the conservation laws. The 
component M of the rotator’s angular momentum L 
along k remains unchanged during a collision, since the 
change in the neutron’s orbital angular momentum is 
perpendicular to k. Hence, with K again designating 
the component of L along the symmetry axis, we 
require the matrix elements for the transitions ( J,K ) 
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—> (J+n, K+k ), which we write as 

M nk = (J+n K+k M\vv(-iK*h)\JKM). (3.21) 

The same system of principal axes as in the low-energy 
calculation will now be used, and we call 0 O the polar 
angle of b with the Z axis or k and 0 O the angle which b 
makes with the symmetry z axis. Recalling that the 
principal xz plane may be chosen to contain the 
nucleus, we have 0 O = tan' 1 (bi/b^). Now consider the 
expansion of the exponential in (3.21) in Legendre 
polynomials P*(cos0o). The addition theorem enables 
us to express each Pi in terms of the polar and azimuthal 
angles of b and k relative to the body axes, namely 
(/3 0 , 0) and (-“0, — ^), respectively. Combining these 
two developments gives 

exp(— iK-b)= £ LM2l+im-i) l ji(Kb) 

1=0 

X £ Y lk (Pofi)U {l ho(o>). (3.22) 

k 1 

The matrix elements (3.21) can be evaluated in the 
Euler angle representation by means of the reduction 
formula (3.17). It is seen that each higher l value in 
(3.22) introduces larger values of \J+n\ in the resolu- 
tion of exp(— iK'b)U (J) KM- Of course, only energetically 
allowed transitions can actually occur. 

The nature of the low-energy approximation is 
revealed through the asymptotic behavior of the 
spherical Bessel functions for small arguments: 

ji(x) « [1 * 3 • • • (21+ 1)]-V[1 - *(2J+3)-W]. 

Since the cross section for a direct scattering collision 
(/= v) involves the square of the magnitude of the 
matrix elements (3.21), the approximation to order 
(Kb) 2 is a consequence of retaining terms up to 1=2 
in (3.22). In fact, this approximation gives precisely 
the cross section (3.18). For example, contributions to 
the elastic scattering come from the / = 0,1,2 and k = 0 
terms with the result that 

(2/+1)-* £ \Mn\*=l-&P+l*btJ- l (J+l')r l K*. 

M^~J 

The inelastic terms in (3.18) come from the 1=1 and 
k= — 1, 0, 1 terms. 

4. INTERFERENCE SCATTERING 

The part of the rotator cross section due to inter- 
ference scattering involving nuclei 1 and 2 is obtained 
by replacing the single expectation value in Eq. (2.1) 
by the pair (X 12 )/+ (X 21 )j. If bi 2 =bi— b 2 , the fi^st 
expectation value is given by 

(X 12 ) j = (2/+ 1 ) -1 Em (/M|exp(zK-bi 2 ) 

X exp[ -it(A+ P 2 +C 2 )] | JM), (4.1) 

where the subscript 2 on the operators B and C indicates 


that they are taken with reference to nucleus 2. In 
the second expectation value the roles of 1 and 2 are 
interchanged. We shall again employ the technique of 
expanding the expectation values in powers of (it) 
and then applying Eq. (3.4). 

The case of a symmetric rotator in the initial state 
| JKM) will be treated to first order accuracy in the 
mass ratios. For taking components we use the system 
of principal axes in which z is the symmetry axis and 
the xz plane is parallel to bi 2 . Hereafter the bar on 
K will be omitted. We let 70 be the angle which 
bi2 makes with the z axis and d=bi 2 . In the ex- 
pressions to come the notation will be simplified if we 
write components in this system as follows: d= (d t ), 
k= (ki) y bi= (u h u 2 ,uz), and b 2 = (v h V 2 ,v d ). In (4.1) the 
quantities c and Q are obtained from (3.10) and (3.11), 
respectively, by inserting for bi the components of 
b 2 . Similarly the tensor R has components Ru=Vj 2 I k ~ l 
+vj?If l y Rij= —ViVjI k -\ Throughout we set 

/i=/ 2 =/. Then from the unitary property of the 
U (J) km we obtain 

(2/+1) -1 Em (JKM | exp (itc ♦ d) | JKM) 

= (1/871^)^ exp(iK*d)dco, 

(2/+1)- 1 Em (JKM\exp(iv:-&)B 2 \JKM) 

(4.2) 

= ( 1/8 ir 2 ) J exp (in ■ d) [ ■ - QuKki - QnKk 2 

m 3 i03 y 

(2/+1)- 1 Em (JKM | exp (iic • d)Ci | JKM) 

(l/S^J exp (iv. ■ d) 2^R mn /c m /c n c^co. 

In (4.2) the integrals are taken over all orientations of 
the rotator. 

The integrals occurring in (4.2) can be expressed as 
a linear combination of integrals which are evaluated 
in the Appendix by the following transformation of 
variables in the integrands. The transformation 
corresponds to a rotation from the original principal 
axes to the system of body axes in which d is parallel 
to the new z axis, i.e., to a rotation through the angle 
70= tan _1 (^i/J 3 ) about the y axis. 

kl~ Ki COS70+K3 siftYO) 

k 2 = k 2 , (4.3) 

KZ= K\ SUI70+/C3 COS70. 

For example, the use of the transformation (4.3) in 
the second integral of (4.2) yields in the notation of 
the Appendix 

(2J+ 1)- 1 Em (J KM \ exp(iK - d)B? \ J KM) 

= (l/S^^-QizK+iici) sinyo 

+hicz cos7ojir(001). 
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After evaluating, the integrals in (4.2) we then have 
(*i 2 )/ to the desired accuracy and it can be expressed in 
terms of jo and j\ alone by means of the relation 
j2(x)=3or l ji(x)— jo(x). The corresponding expression 
for (X 2 i)/ is obtained by interchanging 1 and 2 in the 
previous result. Adding the two expectation values 
and averaging over +K and —K, we obtain finally 
from Eq. (3.4) 

<nnt (1,2) (#)/aia 2 

m d \ h c 2 

= 2 ioM) \ — (biXb 2 ) 2 ioM 

ko 2 dk\ Id? 

+— I" ( *“4 Vb i • b 2 ) (b i X b 2 ) 2 1; i (ad) 

dl\I Ij Id? J 

fcc /1 1 \ r d? 

+“( ) 2u 3 vzji{Kd)+(u 2 2 +V2 2 ) — 

2d\I I z /l d 

x(-i i M)-,i°M)]l . (4.4) 

APPENDIX 

After the angular momentum components are 
allowed to operate upon the basic states in the expecta- 
tion value (3.5), there remains a collection of terms 
each containing the expectation value of a product of 


the form Ki ni K 2 ni K Z ni , where the tti are integers. The k% 
are the components with respect to a set of body axes 
xyz fixed in the rotator of a quantity k directed along 
the Z axis fixed in space. Going over to the representa- 
tion in the Euler angles of the body axes, 

we have 

1 J 

£ {JKMlKf'KfWAJK'M) 

(2/+1) M~j ' 

=— £ f do> UW K M*(a>)UW K 'M 

87T 2 M~J J 

X (Al) 

where d=a>d<pd\[/ sw6d6 and the integration extends 
over all orientations of the body. The sum over M 
can be carried out immediately with the help of the 
unitary property of the matrix U iJ) . The right-hand 
side of Eq. (Al) becomes (I/S^Srk'J (nin 2 n z ) with 

J(tiin 2 n z ) = J Ki ni K 2 n *K Z n *d<j). (A2) 

Expressing the components in terms of the Euler angles 
as in Eq. (3.16) and letting n=fii+n 2 +n z , we can write 
the integral J in the product form 

/ (nin 2 n z ) ~ a n G (nin 2 ) F {n\ w 2 w 3 ) , 

where 


^[i+(-i) ni ][i+(“i) n2 ] r [Kwi+i)]r[§(w 2 +i)] 


and 


G(»i» 2 ) = 27t(— l)" 1 f d\p (cosi^) nl (sin^) n2 = 

r(3(rai+» 2 +2) 

F(tiin 2 n z )= f d$ (cos0) n8 (sin0) ni+n ^ 1 = [l+(— l) n, ]r[J(w 3 +l)]r[J(ni+w 2 +2)J/2r[^(w+3)3. 

J n 


(A3) 


Combining these results gives 


/(«i« 2 « 3 ) = 2^ n | II [l+(~l) n G r [2( w *+l)]J 

1 / 

r Ci(tt+3)]. 

The expectation values encountered in Sec. 4 led to 
integrals of the type 

I(nin 2 n z ) = J* e id * 3 Ki nt K 2 n2 K Z n3 da), 

(A4) 

which can be expressed in the form 


I(nin 2 n z ) =K n G(nin 2 )E(fiin 2 n z ), 

(A5) 

H(nin i n 3 )=j' dp e i ^(l-M 2 ) t " 1 +nil)/ V s 

. (A6) 


-1 


Since G(n\nj) vanishes unless n\ and n 2 are both even, 
we need only consider the case |(wi+w 2 ) equal to an 
integer. When in Eq. (A6) we raise (1— /i 2 ) to an 
integral power, each term contains an even power of 
fx. A typical term, say /x 2m , yields an expression of the 
type 

/ I r2m+n3 ~| 

day £ (A7) 

The sum in (A7) is effectively over only even or only 
odd values of l depending on whether n z is an even or 
odd integer, respectively, because Pi(—fx) = (—1) 1 Pi(jjl). 
The integrals in (A7) are all special cases of a general 
integral formula. 6 The particular integrals (A6) required 
to derive Eq. (4.4) are H(000) = 2j 0 (Kd), J?(001) 
= 2 iji(Kd), F(200) = £T(020) = f[ioM)+i 2 M] and 

g(QQ2)-? [ioM - 2 J 2 M)]. 

6 Bateman Manuscript Project, Higher Transcendental Functions 
(McGraw-Hill Book Company, Inc., New York, 1953), Vol. 1, 
p. 171, Eq. (23). 



